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ABSTRACT 


The Hot Primary Heat Exchanger (HPX), a key component of the 
ThermoAcoustic Life Sciences Refrigerator, consists of a tube and fin design. The 
tubing is bent into a serpentine pattern and overlayed on a screen of copper fins. 
The serpentine pattern results in several flow reversals and complex internal flow 
geometries within the heat exchanger. The fins are not consistently of uniform 
length and generally have heat rejection at both ends. This design results in a 
forced-cooled, single stack cold plate configuration with unequal temperatures at 
each end of the fin. The analysis of this configuration requires a methodology based 
upon the existence of an adiabatic point somewhere along the fin between the prime 
surfaces. Once the location of this adiabatic point is known, the cold plate may be 
treated on the basis of two isolated surfaces having fins with adiabatic tips. The 
goal of this thesis is to provide design analysis and performance predictions for the 
Hot Primary Heat Exchanger (HPX) using numerical analysis of the tube and fin 


arrangement of the HPX. 
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I. INTRODUCTION 


The Hot Primary Heat Exchanger (HPX), a key component of 
the ThermoAcoustic Life Sciences Refrigerator (TALSR), 
consists of a tube and fin design. The tubing is bent into a 
serpentine pattern and overlayed on a screen of copper fins. 
This configuration results in several flow reversals and 
complex internal gas side flow geometries within the heat 
exchanger. The fins are not consistently of uniform length and 
generally have heat rejection at unequal temperatures on each 
end of the fin. This design results in a forced-cooled, 
single stack cold plate configuration. The TALSR uses forced 
convection with an acoustic oscillator providing internal flow 
through channels created within the configuration of the HPX. 

The first formal analyses of the cold plate configuration 
was provided by Mark and Stephenson (1954) and Kraus (1961). 
These analyses provided expressions for the efficiency of the 
cold plate for the case of heat loading on one side. 
Subsequently, Kern and Kraus (1972) looked at the single stack 
cold plate with heat input on only one, and on both sides. 
Incropera and Dewitt (1981) presented a classical textbook 
culmination of extended surface research and general 
conduction analysis for fins of unif.orm cross-sectional area. 
An application of extended surface principles and research was 
used by Garrett (1992) to evaluate the forced convective 
thermal performance of the HPX. Most recently, Pieper and 




Kraus (1995) looked at the cold plate configuration with 
asymmetric heat loading and proposed dividing the plate into 
two fins, each possessing an adiabatic tip for analysis 
purposes. This permits a more accurate representation of the 
performance of the cold plate configuration. 

The goal of this thesis is to utilize the results of 
Pieper and Kraus (1995) to provide design analysis and 
performance predictions for the Hot Primary Heat Exchanger 
(HPX) using a numerical optimization of the serpentine tube 
and fin arrangement in the HPX. 

The serpentine pattern of the copper tubing in the HPX is 
unique when compared to a majority of the tube and fin heat 
exchangers currently in use. A vast majority of current 
applications involve parallel fluid flow in tubes through a 
matrix of fins to provide the required heat transfer. This 
results in fairly uniform base temperatures at each end of the 
fin. In contrast, the serpentine tubing in the HPX results in 
non-uniform base temperatures due to increasing temperatures 
along the length of the tube. Research into the thermal 
performance of a non-uniform base temperature design is 
limited and on-going, however, a successful application of 
this configuration has been achieved on a large scale at the 
WyoDak energy facility in Gillete Wyoming, shown in Figure 1. 

Chapter II of this thesis is used to present fundamental 
heat transfer concepts, and extended surface evaluation 
techniques. Chapter III provides a physical description of 
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ThermoAcoustic Life Sciences Refrigerator (TALSR) and 


.scusses the role of the Hot Primary Heat Exchanger in the 
sat transfer cycle of the TALSR. Chapter IV details the 
imerical simulation techniques used in the performance 
lalysis of the HPX. Chapter V presents a discussion on 
ficiency modeling and the factors affecting extended surface 
■ficiencies. Analysis results and conclusions are presented 
1 Chapters VI, and VII respectively. 







Figure 1: WyoDak Energy Facility in Gillete Wyoming 






II. HEAT TRANSFER FROM EXTENDED SURFACES 


The convective heat transfer rate of a bare surface may 
be increased by increasing the surface area through which the 
convection occurs. This may be accomplished by using surfaces 
that extend from the bare surface into the surrounding fluid. 
Naturally, a considerably larger amount of heat can be 
transferred from or to an extended surface in a given time 
period than from or to a bare surface. The type of extended 
surface most commonly used, is termed a fin. The thermal 
conductivity of the fin material has a strong effect on the 
temperature distribution along the fin and therefore 
influences the degree to which the heat transfer rate is 
enhanced. 

A. BASIC HEAT TRANSFER CONCEPTS 

Conduction is the heat flow mechanism whereby heat is 
transferred by molecular diffusion from one part of a medium 
under the influence of a temperature gradient without a net 
displacement of the particles that compose the medium. It was 
Fourier who proposed that the heat flow is directly 
proportional to the area of the heat flow path and the 
temperature gradient along the path. 


q « A{dT/dx) 


(2.1) 




Insertion of a proportionality constant yields: 

qcond = -kA(dT/dx) (2.2) 

where, 

A = area of the heat flow path 
k = thermal conductivity of the material 
(dT/dx) = change in temperature per unit length 
The minus sign assures a positive heat flow in the presence of 
the required negative temperature gradient. 

Convection is a fluid flow process that results in the 
transfer of heat from or to a confining surface by a flowing 
fluid. The fluid flow may be induced by buoyancy, density 
gradients or through the use of mechanical methods. Those 
methods utilizing mechanical flow generation are termed forced 
convection. 

A second classification of convective heat transfer is as 
either internal or external. In internal flow the fluid is 
constrained on all sides by solid boundaries, as in flow 
through a pipe. In external flow the fluid has at least one 
side extending to infinity without encountering a solid 
surface. Heat flow during convection is directly proportional 
to the temperature difference between the confining surface 
and the surface area over which the process takes place: 

q « A at (2.3) 

Insertion of a proportionality factor yields, 


q^on, = hA AT 


( 2 . 4 ) 



Equation (2.4) is Newton's law of cooling and h is called the 
convection heat transfer coefficient and encompasses all the 
effects that influence the convection mode. 

B. FIN CONCEPTS 

To determine the heat transfer rate associated with a 
fin, the temperature distribution along the fin must first be 
obtained. In the temperature distribution analysis, some 
standard assumptions are made. First, radiation effects are 
neglected. In addition, the fin is assumed to comply with the 
well known Murray (1938) and Gardner (1945) assumptions listed 
in Appendix A. The one dimensional assumption of Murray- 
Gardner is valid because in most extended surface applications 
the fins are relatively very thin compared to their height. 
Thus, the temperature changes in the longitudinal direction 
are much larger than those in the transverse direction and 
the one dimensional assumption is satisfactory. The material 
used in the construction of the fin is characterized by a 
thermal conductivity, k. It is also assumed that the heat 
transfer coefficient due to convection, h, is known. Thus, 
following the standard, steady state fin characterization of 
Incropera and Dewitt (1981), it is found that; 

(k)d(A^(x)dT/dx)/dx - (h) (dA,(x)/dx) (T - T.) =0 (2.5) 

where A^ is the cross-sectional area, and A 3 is the surface 
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area, both of which may vary with x. Equation (2.5) may be 
simplified by defining a temperature excess, 6, as, 

0(x) = T(x) - T, (2.6) 

where because T. is a constant, d0/dx = dT/dx. Substituting 
Equation (2.6) into Equation (2.5), results in; 

()c)d(A,(x)d0(x)/dx)/dx - (h) (dA,(x)/dx)0(x)= 0 (2.7) 

To determine the temperature distribution along the length of 
an individual fin, it is necessary to solve Equation (2.7) for 
the specific fin geometry. 

For the case of an individual rectangular fin of length 
b, as shown in Figure 2, consider that the origin of the 
coordinate is at the fin tip with positive orientation toward 
the fin base. The fin width is a constant, 6, therefore is 
a constant and A, = Px, where A, is the surface area measured 
from the tip to x and P is the fin perimeter. Thus, (dA^/dx) 
= 0 and (dAj/dx) = P = (2L + 26) - 2L and Equation (2.7) 

reduces to; 


d^0/dx^ - m^0 = 0 


( 2 . 8 ) 




Figure 2: Rectangular longituidinal fin of length b. 












= 2h/ka (2.9) 

Equation (2.8) is a linear, homogenous, second-order 

differential equation with constant coefficients. Its general 

solution is of the form; 

0(x) = Cie™ + Cje'™ (2.10) 

By substitution, it is easily verified that Equation (2.10) is 
a solution to Equation (2.8). 

To evaluate the constants and Cj of Equation (2.10), 
it is necessary to specify initial value data; 

0(x = b) = 0b 

and (2.11) 

q(x = b) = qb 

This ma)ces 

0b = Cie"^ + C2e-"*> (2.12) 

and applying qt = )cA(d0b/dx) gives 


qb = k8mL[Cie'^ - Cae'"**] 


(2.13) 



qi, = YoCCje'* - Coe'*] 


(2.14) 


where Y,, = k5inL, is called the characteristic thermal 

admittance of the fin and has the units W/°C. 

It is then a matter of algebra to evaluate the constants 
Cl and Cj such that; 

0(x) = 0bCosh(m[b-x]) - (qt/YJ sinh (m[b-x]) (2.15) 


and 


q(x) = 0bYbSinh(m[b-x]) + qtcosh (m[b-x]) (2.16) 

Now that the temperature excess and heat flow at any point in 
the fin have been evaluated, a convenient method of mapping 
conditions at the fin tip to conditions at the fin base is 
desired. 

For individual fins, Kraus et al. (1978) showed that 
conditions of heat flow and temperature excess (relative to 
the presumed constant and uniform temperature environment) at 
any point on a fin are induced by similar conditions at the 
fin base. This resulted in the development of a linear 
transformation that mapped conditions at the fin tip to 
conditions at the fin base: 
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(2.17) 




where the matrix P is called the inverse thermal transmission 
matrix. Its elements are designated as the inverse thermal 
transmission parameters. A summary of this wor)c is provided 
in Appendix B. 

Applying the development of Kraus et al. (1978) to the 
rectangular fin of Figure 2, it is seen that the inverse 
thermal transmission matrix is given by 


[coshinb (Zg) (sinhmb)l ^2 
Fg(sinhini)) coshrab J 

where = (2h)c6)'^^L, = 1/Yo, and k, L, h, and 6 are the 
thermal conductivity, fin length, heat transfer coefficient, 
and fin thickness respectively. Thus, this matrix can be used 
to map conditions at the fin tip to conditions at the fin base 



costunb 

(siniunb) Fe 1 

. 

rg{sinhmb) 

costuni) 
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In addition to developing the linear transformation 
matrices, Kraus et al. (1978) also proposed that the 
conventional fin efficiency be abandoned and that single fins 
be characterized by a single parameter called the thermal 
transmission ratio, the ratio of the heat entering the fin to 
the temperature excess at the base of the fin. This was later 
called the fin input admittance [Kraus (1982)] and was given 
in the form of a bilinear transformation; 

^Ln = qb/0a = (Ta, + (q^/0a) T 22 ) / (T „ + (qa/0a)Ti2) <2.20) 

The fin input admittance is particularly useful in the 
analysis and evaluation of finned arrays and will play an 
important role in the analysis of the HPX. 

C. FIN ANALYSIS WITH AZISYMMETRIC HEAT LOAOINQ 

Consider the forced cooled cold plate shown in Figure 3 
and observe that the single fin of total height b has been 
subdivided into two fins with fin heights bj and b 2 , to allow 
for the fact that the temperature differences at its opposite 
ends may not be equal. Note that the temperature excess, 0, 
is defined as the difference between the temperature at any 
point on the fin, T, and the temperature of the coolant fluid, 
T., as given in Equation (2.6). 

As shown in Figure 3, the origins of the fin height 
coordinates, Xj and Xj, are taken at the tips of fin 1 and fin 
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2 with a positive orientation from fin tip to fin base. The 
fin heights, b, and b 2 , are chosen such that no heat crosses 
the interface between fin 1 and fin 2, thus creating an 
adiabatic fin tip condition. This is characterized by a linear 
transformation of Equation (2.19) for the longitudinal fin of 
rectangular profile shown, given by; 


M . 

coshmbi Zg(sinhmbi) 



Fg(sinhmb^) coshmbi 



N 

coshmb^ Zg{sinhmb 2 ) 



F;,{sinhmbj) coshinbj 



The cold plate configuration of Figure 3 is subject to the 
continuity and compatibility conditions at the interface 
between fins 1 and 2. These conditions require an adiabatic 
interface between bj and bj and are given in matrix form as; 
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If 0bi ^ 0 b 2 » symmetry will not apply and the fin height s 

hi and bz will not be equal (bj f bj) . 

In order to determine fin lengths b^ and bj, Equations 
( 2 . 21 ) and ( 2 . 22 ) may be expanded using the conditions of 
Equation (2.23). In particular with 0^1 = 0^2 = 0a/ the 

following equations are derived, 

Obi = [cosh mbi]0a (2.24) 

0 b 2 = (cosh mb 2]03 (2.25) 

Equation (2.24) and Equation (2.25) show that the base 

temperature excesses for b, and ba are related through 0 ^, 
thus, a temperature excess ratio (R*) for a given cold plate 
configuration may be defined as: 

R« = 6 bi/ 0 b 2 = [cosh mbj]/[cosh mbj] (2.26) 

The hyperbolic cosines of Equation (2.26) can also be 

represented as exponentials; 

R, = (e"**' + e-'*^)/(e"*’'+ e'"^') (2.27) 

Using the observation from Figure 3, that bj = b - bi, it is 
a matter of algebra to show that; 

R, = (A^ + e^"''‘’=>)/(A[e 2 ”’'“'+ 1 ]) (2.28) 
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A = 


(2.29) 


The value of for a given value of R, can then be found by 
re-arranging Equation (2.28) to provide; 

bj = (l/2in) (ln[A(A - R,)/(ReA - 1)]) (2.30) 

An evaluation of using Equation (2.30) confirms that if R,=l 
because 0^1 = 0 b 2 , then because of symmetry, must equal half 
of the total fin height. In addition , the value of R, for 
b2=0 can be found directly from Equation (2.28), 

(R«)™ax = (A" + 1)/(2A) (2.31) 

and similarly for b 2 =b, 

= (2A)/(A2 + 1) (2.32) 

If R, is not within the domain of values specified by 
Equations (2.31) and (2.32) the fin is then treated as a 
single fin without an adiabatic point. However, for the HPX 
model, the Rj values are always within this domain and the 
single fin analysis is not required. 

Once bj is known, bi is easily found through a simple 
subtraction procedure (bj = b - b 2 ) . With both bi and bj 
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known, the cold plate may be treated on the basis of two 
isolated surfaces having fins with adiabatic tips, treating 
each fin individually as though the other were not present. 
One surface is governed by bj with 0bi specified, and the other 
is governed by bj with specified. Thus, there are two 
entities, each with a prime surface and each with a fin, and 
these surfaces may be treated individually as if the other 
were not present. 

The performance of each prime surface, and fin (bi or bj) 
combination depends on the total input admittance of the pair. 
The total input admittance, Yi„, is just the sum of the prime 
surface and fin input admittances, 

Vin = (2.33) 

where the subscripts, p and f refer to the prime surface and 
the fin, respectively, and where both prime and fin base 
surfaces are operating at Ot,. 

The prime surface input admittance is determined by its 
convective dissipation 


Yin,p = hSp 


(2.34) 


where as seen from Figure 4, 


Sp 


(W - 6)L 


(2.35) 
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which accounts for the footprint occupied by the fin. 

The fin input admittance is derived from the fin p-matrix 
given by Equation (2.18). Because the height b^ is based upon 
the determination of an adiabatic point, is established 

by realizing that for qa = 0 Equation (2.20) becomes; 

^in = qb/0« = ^ 21/^11 = YoSinh(mb)/cosh(mb) (2.36) 

which for the fin governed by bj, becomes 


Yin,f = 'io tanh(mb2) (2.37) 

Then, if the temperature excesses are specified so that 
leads to the establishment of bj, the heat dissipation is 
obtained from; 


qb = l^tb9b (2.38) 

This section provides the basis for an accurate thermal 
performance model needed in the simulation of the heat 
transfer process in the TALSR discussed in the next chapter. 
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III. 


THERMOACOnSTIC LIFE SCIENCES REFRIGERATOR 


The ThermoAcoustic Life Sciences Refrigerator (TALSR), 
developed at Naval Postgraduate School, was motivated by the 
desire to replace the Freon 512 vapor compression 
refrigeration system currently used on board the Space 
Shuttle. The TALSR provided a safer alternative due to the 
fact that it does not use chlorof luorocarbons (CFCs) and 
therefore could not potentially contaminate the small confined 
area of the Space Shuttle. In addition, the TALSR also 
provides the potential for higher reliability, over the 
current system, because it has no sliding seals and thus 
reguires no lubrication. The ThermoAcoustic Life Sciences 
Refrigerator uses a complex thermal transport subsystem to 
remove heat from the Insulated Refrigeration Enclosure (IRE) 
of the Space Shuttle. 

A. THERMAL SYSTEMS OVERVIEW 

Heat is removed from the IRE by a Liguid-Air Heat 
Exchanger fabricated by the Modine Manufacturing Company. It 
is then transported through two thermoacoustic heat pumps, 
connected in series, by a mechanically pumped cold heat 
exchange fluid. The internal architecture of each heat pump 
consists of a Cold Primary Heat Exchanger (CPX), a Hot Primary 
Heat Exchanger (HPX), and an Electrodynamic Driver. The 
Electrodynamic Driver acoustically oscillates an internal 
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working fluid in each heat pump thereby creating internal gas 
side flow for both the HPX and the CPX. It is known that the 
heat transfer between surfaces can be enhanced using 
longitudinal acoustic waves [Vainshtein (1995)]. This internal 
gas side flow is used by the Cold Primary Heat Exchanger to 
remove heat from the cold-side exchange fluid, and transport 
it to the HPX. The HPX then transfers the heat removed from 
the cold-side exchange fluid into a secondary cooling loop. 
This secondary cooling loop transports the heat to the 
internal environment of the Space Shuttle via an Experimental 
Heat Exchanger (EHX). A schematic diagram of the series fluid 
flow through the TALSR thermal transport subsystem is shown in 
Figure 5 [Garrett (1992)]. 

B. THE HOT PRIMABY HEAT EXCHANGER (HPX) 

The Hot Primary Heat Exchanger uses an acoustically 
oscillating medium to transport heat and must therefore be 
analyzed differently from most compact heat exchangers. The 
most significant difference between standard compact heat 
exchanger analysis and acoustical heat exchanger analysis is 
the fact that the acoustically oscillating gas parcels only 
move a limited distance before reversing their direction of 
flow. The consequence of this periodic flow reversal is that 
an increase in the effective surface area for heat transfer 
cannot arbitrarily be achieved by simply increasing the length 
of the heat exchange surfaces (fins) in the direction of flow. 
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Garrett (1992) has shown that the effective area over which 
heat transfer will take place is limited by the peak-to-peak 
excursions of the gas parcels over the heat exchange surface. 
Thus, heat exchange surfaces, which have a length in the flow 
direction equal to the peak-to-peak displacement of the gas, 
provide the maximum effective surface area for heat transfer. 

The HPX utilizes copper tubing bent into a serpentine 
pattern and then soldered or furnace brazed on a screen of 
copper fins, as shown in Figure 6 [Garrett (1992)]. In this 
configuration, the fins are not of uniform length and for a 
majority of the fins, there is heat rejection at both ends. 
This design results in two parallel paths by which the fluid 
contained within the copper tubing can remove heat from the 
acoustically oscillating gas. The primary path is through the 
fins which are bonded to the tube. The fins represent the 
primary path because they have a greater surface area than the 
outside diameter of the tubing and because they are designed 
to have a length in the flow direction approximately equal to 
the peak-to-peak displacement of the gas. The secondary path 
is due to the direct convection of the gas around the tubing. 
This path also provides some useful heat transfer and requires 
consideration in the thermal analysis of the HPX. A schematic 
representation of the thermal resistances in the primary heat 
exchanger energy flow path is shown in Figure 7 [Garrett 
(1992)]. The tubing for the HPX is standard, circular cross- 
section, soft copper tubing, with an outer diameter (OD) of 
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6: HPX design configuration [Garrett (1992)] 
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Tube 


Figure 7: Schematic representation of thermal resistances in the 
HPX flow path. [Garrett (1992)] 
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0.635 cm and a wall thickness of 0.076 cm. The HPX has an 
inner diameter of 11 cm. Table 1 provides a list of the 
geometrical constraints and heat transfer characteristics used 
in the analysis of the HPX. 

The physical design of the HPX can therefore be modeled 
as a forced-cooled, single stack cold plate configuration with 
unequal temperatures at each end of the fin. Thus, an 
analysis of the HPX can be completed using the derived 
relationships for cold plates with axisymmetric heat loading. 


Symbol 

Description 

Value 

Reference 

Cp 

specific heat capacity 
of water 

4.182 

kJ/kg-K 

Incropera and 
Dewitt (1981) 

5 

Fin thickness 

0.0152 cm 

Garrett (1992) 

hf. 

Fluid heat transfer 
coefficient 

12600 W/m^-K 

Garrett (1992) 

k 

thermal conductivity of 
copper 

401 W/m-K 

Garrett (1992) 

L 

Fin width (channel 
length) 

0.3175 cm 

Garrett (1992) 

Pr 

Prandtl number 

0.7068 

Garrett (1992) 

Re 

Reynolds number 

1900 

Garrett (1992) 

Rinner 

Inner radius of copper 
tubing 

0.2413 cm 

Garrett (1992) 

R=ut,r 

Outer radius of copper 
tubing 

0.3175 cm 

Garrett (1992) 


mean dynamic viscosity 

190 X 10'^ 
N-s/m^ 

Incropera and 
Dewitt (1981) 

Ms 

surface dynamic 
viscosity 

184 X 10'^ 
N-s/m^ 

Incropera and 
Dewitt (1981) 


Table 1: A summary of the physical constraints and heat 
transfer properties used in the analysis of the HPX 
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C. COMPUTER MODELING OP THE HEAT TRANSFER PROCESS 

Once an accurate thermal model has been obtained, an 
efficient and reliable method of analyzing the heat transfer 
processes occurring within the heat exchanger is reguired. 
This is accomplished through the use of Steady State Thermal 
Analyzer software version 2.2 provided by Intercept Software. 
The model builder provides the user with the means to model 
the physical configuration of interest in order to provide the 
thermal analyzer with an input file. The thermal analyzer 
then takes the input file and produces an output file 
containing a summary of the temperatures within the 
configuration. 

The model building process begins with a drawing of the 
configuration and a subdivision of it into small but finite 
subvolumes. Each subvolume is presumed isothermal and the 
centers of each subvolume are then representative of the 
entire subvolume. These centers are referred to as nodes and 
are connected to adjacent nodes through branches consisting of 
various forms of thermal conductance. These conductance forms 
are dependent on the mode of heat transfer between adjoining 
pairs of nodes. The various heat transfer modes available to 
the user include conduction, laminar free convection, 
radiation, forced convection, and fluid flow. 

Each node can also be connected to a constant 
temperature, a constant heat input, or a temperature dependent 
heat input with an appropriate tag. Once a sketch 



representation of the nodalized model is complete, the 
comprehensive node connection data is ready to be input into 
the model builder program. The input for the model builder 
program begins with node 1, and the user is asked for 
information for each connecting node with a greater number. 
For each connection to an adjacent node, the user must specify 
the mode of heat transfer with a tag number. After specifying 
the connecting node and tag, the user is queried as to the 
whether the conductance will be calculated, input directly, or 
is the same as an earlier branch. Conductance values are 
required for each node connection to any higher number node. 

Subsequent use of the thermal analysis software results 
in the writing of n node equations in n unknown temperatures 
where the nodes are connected by the appropriate thermal 
conductances. The general solution strategy for these n 
equations is to then use the nodal conductances of the model 
to form a set of heat balance equations. In matrix form, this 
set of heat balance equations has the general foi-m; 

[K] [T] = [B] (3.1) 

where K is the matrix of conductances, T is the node 
temperatures, and B is a matrix composed of constant 
temperature heat sinks, and/or heat inputs. Many of the 
conductances are linearized forms of nonlinear expressions for 
heat transfer by natural convection and other similar heat 
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transfer inodes, therefore, the various terms in the K matrix 
may themselves be functions of temperature. Solution of the 
T matrix therefore, must be by iteration. Node temperatures 
obtained after each iteration are used to update the 
temperature dependent terms in the X matrix. The thermal 
analyzer uses a Cholesky factorization [Hamming (1973)] of the 
K matrix to perform this iteration. This iterative solution 
continues until the change in nodal temperatures between 
successive iterations is smaller than a user-specified error 
criteria. Once the iterative solution is obtained, the thermal 
analyzer writes the temperatures to an output file, where they 
can be read and analyzed by the user. 



IV. PERFORMANCE ANALYSIS 


The thermal analysis of the TALSR requires an adequate 
model capable of determining both the heat transfer 
characteristics for the complex flow geometries within the 
HPX, and the values of temperature along the serpentine 
pattern of the copper tubing. The thermal conductivity of 
copper, k, the fin thickness, 6, and the channel length, L, 
are all easily determined from the given material and 
geometrical considerations of the HPX. Therefore, to complete 
the thermal analysis, the convection coefficient, h, is 
required to determine the total input admittance, and the 
conductance matrix is needed to map base temperature excesses 
along the length of copper tubing. These remaining values are 
determined using numerical analysis techniques. 

A. HOT PRIMARY GAS SIDE HEAT TRANSFER COEFFICIENT 

DETERMINATION 

The convection heat transfer coefficient, h, encompasses 
all the effects that influence the convection mode. It 
depends on conditions in the boundary layer, which are 
influenced by surface geometry, the nature of fluid motion and 
many of the fluid thermodynamic and transport properties. By 
considering the magnitude of the factors affecting the heat 
transfer coefficient, h, an appreciation for the complexity in 
determining its value is obtained. In simple flow situations, 
solutions for h are readily effected mathematically, however. 
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for situations of complex flow geometries, such as those in 
the TALSR, the more practical approach involves calculating h 
from empirical eguations. The particular form of these 
eguations is obtained by correlating measured convection heat 
and mass transfer results in terms of appropriate 
dimensionless groups. The development of the dimensionless 
group used to determine h was performed by Incropera and 
Dewitt (1981), and resulted in the derivation of the Nusselt 
number (Nu). The dimensionless Nusselt number provides a 
measure of the convective heat transfer occurring at the 
surface. For a prescribed geometry, the Nusselt number is a 
universal function of x‘, the Reynolds number, and the Prandtl 
number, where x" is the dimensionless length of the fluid 
along the channel. A detailed analysis of these factors will 
determine the exact correlation to be used in determining the 
Nusselt number. For the given flow geometry and 
considerations of the HPX, the Reynolds number is given by 
Garrett (1992) to be 1900, therefore, flow through the coolant 
channel remains in the laminar region. In the absence of 
strict correlations for oscillatory flow convection 
coefficients, a steady flow Nusselt correlation was assumed to 
provide a suitable representation of heat transfer in the 
channel. Sieder and Tate (1936) proposed a suitable 
correlation for laminar flow in tubes and ducts 
([ (ReoPrD/L)^^^) (^/M,)‘' *‘] ^ 2) given by; 
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.86 (RePr/(L/D))''^(M//ij““ 


(4.1) 


NU = 1 

where ^ and /i, are the average dynamic viscosity and the 
surface dynamic viscosity respectively. From knowledge of the 
Nusselt number, the convection coefficient may now be found 
using; 

Nu = (hdg/kp) (4.2) 

where d^ is the effective diameter of a non-circular duct 
through which the fluid passes. Replacing D with d, in 
Equation (4.1) and combining with Equation (4.2) results in 

h= (Nu(kp)/dJ = (l.86kp/dj (RePr/(L/de))’'^(M/Ms)°-” ( 4 . 3 ) 

Then using the values of L, Re and Pr specified in Table 1, 
and realizing that for air (^i/^i,)® ” - 1 over the temperature 
range of interest for the HPX it is easily shown that; 

h = (7.9 X 10-^) (4.4) 

where h is given in Watts/cm^-K. The effective diameter, d,, 
is defined as; 

d, = (4A,/P) (4.5) 
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where A,, and P are the flow cross-sectional area and the 


wetted perimeter respectively. A rectangular duct flow 
geometry is assumed for the TALSR based on the large radius of 
curvature of the HPX and the small spacing between fins, thus 
as shown in Figure 4; 


d, = 22b/(b + 2) ( 4 . 6 ) 
where 2 = (W - 6)/2, and b = bi + bj. The value of b varies by 
location within the HPX, as shown in Figure 5, and results in 
corresponding changes in the convection coefficient for each 
location. Therefore the value of h in Equation (4.4) must be 
considered independently for each location along the copper 
tubing. This significantly increases the complexity of the 
analysis. 

B. TEMPERATURE EXCESS DBTERNINATIOM 

Determination of temperature excess variation with 
location along the length of copper tubing is also required to 
adequately complete the analysis of the TALSR. An accurate 
model of the heat transfer processes shown in Figure 7 is 
required to determine this variation. 

Application of the model builder program THANSS, located 
within the Steady State Thermal Analy 2 er software, to the 
TALSR is accomplished by dividing the length of the copper 
tubing surface into 144 subvolumes of varying si 2 e based on 
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location. Additionally, another 144 nodes are assigned to 
corresponding fluid nodes within the pipe to account for 
variations in fluid temperature as it flows through the 
tubing. A summary of nodal assignments based on location is 
given in Appendix C. An additional two nodes are required to 
complete the model. The first node, a constant heat source 
node, is assigned to the air flowing through the coolant 
channels created by the tube and fin design. The temperature 
of the air is considered uniform and constant at 40°C. The 
second node is used to represent the liquid input temperature 
at the inlet to the HPX. This temperature is also assumed 
uniform and constant at a temperature of 25°C. 

The general form of the heat flow equations used for all 
modes of heat transfer considered within the HPX is; 

q = K(Ti - Tj) ( 4 . 7 ) 

where K, the conductance, varies in form depending on the mode 
of heat transfer considered. For the HPX model, four forms of 
thermal conductance are required to define the various modes 
of heat transfer occurring within the HPX. 

The first conductance is represented as Ki and is given 

by; 


K, = KA^^/aL 


( 4 . 8 ) 





where k is the thermal conductivity, is the cross sectional 
area, and aL is the length between adjacent nodes. This 
conductance is used to define heat flow between adjacent nodes 
on the surface of the copper tubing, therefore, the cross 
sectional area available for heat transfer, A^, is given by: 

A^ = n (Router - R\nn,r) = 0.134 Cm^ (4.9) 

The thermal conductivity of copper (k = 401 W/m-K) was used 
for all calculations involving Kj. 

The second conductance is used in the calculation of the 
convective heat transfer between nodes on the surface of the 
copper tubing and corresponding liquid nodes located within 
the copper tubing. It is denoted as and is given as: 

K2 = li,A, (4.10) 

where h„ is the convective heat transfer coefficient due to 
water flow through the tube, and A, is the surface area over 
which the heat transfer occurs. The given value of hp for 
turbulent water flow through the tube is given by Garrett 
(1992) as 12600 W/m^-K. The surface area over which the heat 
transfer occurs is based on the inner radius of the copper 
tubing and is given by: 

A, = 2nRi„„erAL (4.11) 


36 



where aL 


identical 


that used in determining K,. 


The third conductance form, K 3 , is used in the 
determination of the convective heat transfer between adjacent 
fluid nodes within the tube. It is given as: 

K3 = WeCf (4.12) 

where Wp is the fluid weight flow and c^ is its specific heat 
capacity. This value was calculated for a variety of flow 
rates (w^ = 68 kg/hr, 149.7 kg/hr, 272.2 kg/hr, 362.9 kg/hr) 
through the tubes, and the thermal performance of each flow 
rate was analyzed. A constant value of 4.182 kJ/kg-K was used 
for the specific heat capacity of water in all calculations of 

K3. 

The fourth conductance term necessary to analyze the 
thermal model is used in the calculation of the heat transfer 
from the surface of the tubing, to the surrounding environment 
via both the prime and fin surfaces. This term is the most 
difficult of the conductance forms to determine, and requires 
a return to the principles presented in the single stack cold 
plate analysis. Equation (2.33) is written in the form of the 
general heat equation used by the thermal analysis software in 
its calculations where; 
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( 4 . 13 ) 


= Y,„ = Yi„,p + Yi„,£ 

0a = (Ta - Ta) = (Ti - Tj) (4.14) 

Therefore, the fourth conduction term is the sum of the prime 
surface and fin input admittances, and the temperature 
difference is given in terms of a temperature excess at the 
base of the fin. Applying this conductance form to the HPX 
requires consideration of the length, aL, between adjacent 
nodes. The length between adjacent nodes determines both the 
number of fins, and the amount of prime surface associated 
with each node. The aL in the analysis model is varied with 
location along the copper tubing as shown in Appendix C. Once 
aL is determined, the resulting number of fins, n, and prime 
surface area may be used with the determined value of h to 
give; 

K, = hn(W - 6)L + nYptanh(mb 2 ) (4.15) 

where n = aL/W. 

Once the model provides the required data to satisfy the 
K-value equations, the solution of the thermal conductance 
matrix is straightforward and ideally suited to a spreadsheet 



application. The calculated K-values are then used as input 
to the thermal analyzer to form a series of node connection 
data. The spreadsheet calculations and resulting node 
connection data for each case considered is provided in 
Appendix D. In the HPX model, the tube water inlet 
temperature and the oscillating air temperature are assigned 
constant values. Additionally, as mentioned earlier, 144 
nodes were assigned to locations along the copper tubing 
(nodes 1-144) and 144 corresponding nodes (nodes 145-288) were 
assigned to the centerline liquid temperature within the 
tubing. 

In the initial development, it was assumed that all fin 
pairs were shared equally (0bi = 0b2) • This condition is the 
easiest to model, because fin sharing node pairs have 
equivalent values of K, assigned to them. This model results 
in the assigning of 46 distinct K-values and is used in the 
establishment of the initial temperature excess to be used in 
subsequent iterations. The computed temperature excess 
variation with location was then used in the calculation of 
the nodal fin sharing pair R* values using Equation (2.26). 
Once the values were known, the values of bj and bj were 
calculated using Equation (2.30). The assignment of specific, 
non uniform fin lengths ( b^, bj) along the tubing surface, 
required the value of K, to be modified for each node. This 
was required based on the fact that now unique values of 
existed for each node. Applying the K-value equations to each 
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node, for this model, resulted in the generation of 252 unique 
values of K,. These values were then entered into the thermal 
analysis program. The thermal performance for a variety of 
flow rates (Wp = 68 kg/hr, 149.7 kg/hr, 272.2 kg/hr, 362.9 
kg/hr) was also analyzed for this new model. These analyses 
provided adequate data for a performance prediction of the 
TALSR. 
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V. EFFECTS OF TEMPERATURE EXCESS RATIO VARIATIONS 


The base temperature excess (0b) , as defined earlier, is 
the difference between the temperature at the base of the fin 
and the temperature of the coolant fluid, given by; 

0b = Tb - T. (5.1) 

It represents the maximum driving potential for convective 
heat transfer across the fin. Hence, the maximum rate at 
which a fin could dissipate energy is the rate that would 
exist if the entire fin surface were at the base temperature. 
However, because any fin is characterized by a finite 
conduction resistance, a temperature gradient will exist along 
the fin and the rate of energy dissipation will be less than 
the maximum rate. Thus, a convenient measure for the thermal 
performance of a particular fin is provided by the fin 
efficiency, t), defined as [Incropera and DeWitt (1981)]; 

1 = qf/qideai (5.2) 


where 


qf - qb = ^o0b teM:ih(mb2) 
qideal = hPb0b 


(5.3) 


for a rectangular fin of uniform cross sectional area with an 
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adiabatic tip. The equation for qideai represents the heat 


convectively transferred from a fin of height b to the 
environment if 0 = 0t, at all points along the fin. Combining 
Equations (5.2) and Equation (5.3) results in a conventional 
fin efficiency given as; 

t] = tanh(mb)/mb (5.4) 

where m is as given in Equation (2.9), and b is the total 
height of the fin. The application of this parameter to cold 
plates with asymmetric heat loading, such as those modeled in 
the HPX, will result in individual efficiencies being 
identified for each fin subdivision (bi and bj) . In addition, 
an overall efficiency for the total fin height b is also 
calculable. 

A simplified cold plate configuration will be used as the 
model for this analysis in order to provide clarity to the 
development. The model will consist of a single fin with a 
total height b, and )cnown heat transfer characteristics for 
the surrounding flow geometry. Therefore, the only required 
data to complete the thermal analysis of the model is the 
prime surface temperature on either end of the fin. Once this 
value is known, R,, bi and bj are easily determined using 
Equations (2.26) through (2.28). 

A conventional derivation of the efficiencies for each 
fin subdivision is easily determined by substituting the 
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values of bj and bj in for the total fin height, b, in Equation 
(5.3). A conventional overall efficiency of the entire fin 
height, b, can also be calculated using 

»love.all = (qactaalM + qactualb2 ) / ( qidealbl + qidealbj) (5.5) 

which algebraically reduces to 

novabau = ^bl/d + (b2/b,R«)) + T|,2/(l + (biR^/b^)) (5.6) 

The value of R« used in determining Equation (5.6) is based 
upon the average temperature change across the fin as given by 

Rj = (aT/2 + - T,)/(-aT/2 + - TJ (5.7) 


at = (T,i - T,2) 

and T^^, = (T^i + T^z)/2 (5.8) 

The use of this definition for R«, contrary to one based 
solely on Tbi or Tta, provides an accurate representation of 
the response in fin efficiency due to variations in base 
temperature excess across the fin. 

These variations in fin efficiency were then plotted 
versus the difference in base temperature excesses as shown in 
Figure 8. 

A more practical definition for fin efficiency would 
result in the efficiency of each segment increasing with the 
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Conventional Efficiencies vs. Temperature Excess Difference 



'igure 8: Graphical representation of conventional efficiency variatic 
lifference_ 


changing 


I temperature 








length of the respective fin heights. For this development th< 


overall efficiency is defined as 

tioverall = + 12 (5.9) 

where t], and TI 2 are given by 

’ll = qa/qidaai = Yo^bi tanh(mbj/qideal (5.10) 

and i = 1 or 2. The value of qideai for this development is 
given as 

qidaal = hPb0bi (5.11) 


where 0^1 is chosen as the maximum value of the magnitude of 
either 0 m or 0^2• For example, for R, 2 1 it follows from 
Equation (5.11) that qidaai = hPb0M. After substitution of this 
result into Equation (5.10) and combining with Equation 
(5.11), it is found that for R, 2 1 overall efficiency is 
given as 


(tanh(mbi) / (mb)) + (tanh(mb2) / (mbR,) ) (5.12) 





similarly for i 1 


noveraii = (Re(tanh(mbi)) / (mb) ) + (tanh(mb 2 ) / (mb) ) (5.13) 

These equations are seen to be equal for the case where Rj = 
1. A graphical representation of the relationship between the 
individual fin efficiencies and the practical overall 
efficiency is presented in Figure 9. 



I Practical Efficiencies vs. Temperature Excess Difference 
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Figure 9 ; Relationship of practical efficiencies with variations in base temperature excess 
difference 
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VI. RESULTS 


An accurate design and performance analysis of the HPX 
requires critical assessment of the thermal analyzer output 
data, and representation of that data in a format that is 
easily understood and interpreted. 

A. PROGRAM OUTPUT 

The thermal analyzer output data for the simplified model 
provides baseline temperature excess variations along the 
surface of the copper tubing and liquid centerline 
temperatures. The temperature output, for the nodes assigned 
to the tubing surface (1-144), indicated variations in 
temperature along the length of the tubing. This is due to 
variations in available heat transfer area, based on the 
assignment of node locations and resultant fin allocation. 
Additionally, the liquid centerline temperature showed a 
gradually increasing trend from inlet to outlet which is 
consistent with the expected results. 

Analysis of the complex model resulted in similar 
temperature variations for both the tubing and liquid 
centerline temperatures, with only slight differences in 
actual temperature values from those obtained using the 
simplified model. Complex model analysis for different flow 
rates indicated that temperature values along the surface of 
the tube, and throughout the liquid in the tube were largely 
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dependent on the mass flow rate. For an increasing mass flow 
rate, this dependence resulted in lower tube surface 
temperatures, fluid outlet temperature, and overall change in 
fluid temperature as it passed through the heat exchanger. A 
summary of the output nodal temperatures for each model 
considered is found in Appendix E. The results were consistent 
with expectations and are presented graphically in Figures 
(10) and (11). In addition, the thermal analysis provided the 
necessary data to calculate the overall heat transfer for each 
mass flow rate, using; 


q = (din/dt)CF aT (6 .1) 

where Cf is the specific heat capacity of water (4.182 kJ/kg- 
K). The heat transfer rate is plotted for each mass flow rate 
and is presented in Figure (12). 
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Tube surface temperature variation (Celsius) 








Coolant Temperature Difference (Celsius) 


Coolant Temperature Difference vs. Mass flow rate 



Figure 11: Graphical representation of Coolant Temperature 
difference variation with mass flow rate. 
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Heat Transfer vs. Mass flow rate 
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VII. SOMMARY AND CONCLUSIONS 


The results of this investigation into the thermal 
performance of the HPX leads to several conclusions. These 
conclusions are broken down into the following topics: design 
effectiveness, effectiveness of applying single stack cold 
plate analysis, effectiveness of computer modeling, and 
performance analysis validity. 

A. DESIGN EFFECTIVENESS 

Design effectiveness consideration was based solely on 
the physical configuration of the HPX. The use of copper 
tubing bent into a serpentine pattern and then soldered or 
furnace brazed on a screen of copper fins results in a small 
and inefficient contact surface between the tubing and the 
fins. A better design and the one actually used to model the 
HPX would inlay the copper tubing into the screen of copper 
fins. This would result in a much higher effective surface 
area for heat transfer than the current configuration. The 
success of applying a serpentine tubing pattern for this 
application as opposed to the more conventional parallel 
pattern was not analyzed and provides a basis for future 


thesis work. 





B. EFFECTIVENESS OF SINGLE STACK COLD PLATE ANALYSIS 

An analysis based on the work of Pieper and Kraus (1995) 
indicated distinct fin sharing characteristics along the 
length of copper tubing based on the temperature excess ratio 
across each shared fin. This resulted in small variations in 
computed temperatures along the length of the tubing compared 
to a model subjected to equal operating temperatures on the 
right and left base surfaces. The magnitude of this variation 
is dependent on the scale of the application and could result 
in large values for a large scale application similar to the 
WyoDak facility in Gillette, Wyoming. Thus, a single stack 
cold plate analysis would be very useful if the analysis 
required very accurate representations of temperature along 
the length of copper tubing in a large scale application. The 
variation in the magnitude of the temperature variation based 
on the scale of the application provides a basis for further 
research. 

C. EFFECTIVENESS OF COMPUTER MODELING 

Computer modeling provided a convenient and reliable 
method of calculating the multiple equations created as the 
result of a detailed analysis of the HPX heat exchanger. The 
use of computer modeling provides flexibility in the analysis 
of various heat transfer parameters. The model used for the 
analysis of the HPX consisted of air flowing over fins. 
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heating water flowing through the copper tubing. Once these 
initial inputs are set up and the requisite files built within 
the thermal analysis software, it becomes simply a matter of 
replacing the initial model values with new values to evaluate 
a new system. Therefore, various working fluids with widely 
varying properties and different internal heat exchanger 
geometries can be analyzed very easily. This permits the 
convenient application of specific conditions for a given 
application, thus saving money and laborious hours in the lab 
building and evaluating various configurations for a given 
application. However, once a given application is chosen 
through the use of computer modeling, extensive testing and 
empirical data collection is required to validate the computer 
generated results. For the HPX model analyzed in this thesis, 
the collection of empirical results through testing and 
subsequent comparison with the computer generated results is 
still required. 

D. VALIDITY OF FERFORHANCE ANALYSIS 

The performance analysis of the HPX required the 
determination of the convection coefficient, h, and the 
conductance matrix. These values were determined using 
numerical analysis techniques. 

In the absence of strict correlations for oscillatory 
flow convection coefficients, a steady flow Nusselt 
correlation was assumed to provide a suitable representation 





of heat transfer in the channel. This is a major assumption, 
but not without merit based on the argument provided in the 
thermal performance section. However, a more representative 
analysis could be made through the application of a Nusselt 
correlation based solely on acoustically oscillated flow 
analysis. The derivation of these correlations is currently 
being pursued through a variety of research, including work 
conducted at Naval Postgraduate School. 

The conductance matrix was used to determine the 
temperature distribution along the length of the copper 
tubing. This distribution of temperatures was reasonable and 
behaved predictably with changes in the mass flow rate of 
water through the tubing. These temperatures were then used 
to calculate values for the temperature excess ratio, and 
subsequently determine the respective lengths of bj and bj. 
A new temperature distribution was then found based on this 
cold plate analysis. This provided an adequate representation 
of the temperature variation along the length of the tubing, 
however a more accurate distribution could be acquired by 
repetition of the analysis for each new set of temperature 
distribution data acquired. This could continue until a user- 
specified variation limit between successive attempts was met, 
thus satisfying a predetermined level of accuracy. 
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APPENDIX A. MURRAY-GARDNER ASSUMPTIONS 


The Murray-Gardner assumptions are: 

1. The heat flow in the fin and its temperatures remain 
constant with time. 

2. The fin material is homogenous, its thermal conductivity is 
the same in all directions, and it remains constant. 

3. The heat transfer coefficient to the fin is constant and 
uniform over the entire surface of the fin. 

4. The temperature of the medium surrounding the fin is 
uniform. 

5. The fin thickness is so small compared with its height that 
temperature gradients across the fin thickness may be 
neglected. 

6. The temperature at the base of the fin is uniform. 

7. There is no contact resistance where the base of the fin 
joins the prime surface. 

8. There are no heat sources within the fin itself. 

9. The heat transferred through the outermost edge of the fin 
is negligible compared with that leaving the fin through its 
lateral surface. 

10. Heat transfer to or from the fin is proportional to the 
temperature excess between the fin and the surrounding medium. 
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APPENDIX B. TRANSFORMATION MATRIX DEVELOPMENT 

A general application of linear, homogenous, second-order 
differential equation theory dictates that Equation (2.7) 
posses two independent solutions. These solutions, and the 
subsequent development of a linear transformation matrix to 
map conditions at the fin tip to conditions at the fin base 
were performed by Kraus et al. (1978) . The solutions, 
designated A.i(x) and A, 2 (x), must satisfy the initial conditions 
at the base of the fin where x = b. 

li(b) = l; Xi’ (b) = 0, 

X2(b) = O; X2'(b) = (l/kA,(b)) (B.l) 

where the prime indicates a first derivative. The heat flow, 
q(x), is always ta)ten as positive from base to tip. 
Therefore, for longitudinal fins, q(x) is given by; 

q(x) = kA^(x)(de(x)/dx) (B.2) 

Thus, the solutions and Xj can be used to assemble the 
expressions for the temperature excess 0(x) and heat flow 
q(x), at any point in the fin in terms of 0b and q^ at the fin 


base; 





0(x) = 0bXi(x) + 


(B.3) 


q(x) = kA,(x) [0,A,'(X) +qt,^ 2 '(x)] (B.4) 

In matrix form, Equations (B.3) and (B.4) become; 


ElJil 


[1 0 1 


^2 (X) ■ 

[0 

(x) 

X2(x) 


(B.5) 


The thermal transmission matrix is the linear transformation 
generated when x is set equal to a, where a equals either the 
fin height or zero depending on the origin of the height 
coordinate; 



(B.6) 


where 

Yii = ^i(a) 

Yu = ^2 (a)' 

= )cA,(a)li'{a) 

Y 22 = kA,(a)l 2 '(a) (B.7) 
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The elements of the thermal transmission matrix are called the 
thermal transmission parameters. To represent conditions at 
the fin base in terms of conditions at the fin tip, it is seen 
that ; 



where the matrix p is the inverse of the thermal transmission 
matrix and is called the inverse thermal transmission matrix. 
Its elements are designated as the inverse thermal 
transmission parameters. 


63 





6A 



APPENDIX C. NODAL MODELING DATA 



















APPENDIX D. SPREADSHEET K-VALDE DETERMINATION AND NODE 
CONNECTION DATA 


The general procedure for each case is as follows: 

1. Use node location and applicable equations to determine fin 
length associated with that particular node. 

2. Determine K-Values using spreadsheet application. 

3. Input K-Values into user friendly spreadsheet that associates 
each K-Value with a particular node and branch. 

4. Once K-Values are input into thermal analysis software, an 
easily read summary of the individual branches and their 
associated conductances is produced. 
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TASS GENERAL INPUT MENU - SI Units 


(1) Case Title: 

TALSR(METRIC)- RUN 1. SIMPLE MODEL CASE 

(2) Nodes 288 

(3) Constant Temperatures 2 

(4) Unique Exponents 0 

(5) Temperature Dependent Conductances 0 

(6) Temperature Dependent Heat Inputs 0 

(7) Computational Accuracy .0100 

(8) Starting Temperature 25.0 

Are these inputs correct (Y/N) ? Y 
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TASS GENERAL INPUT MENU - SI Units 
(1) Case Title: 

TALSR(METRIC)-RUN 2. COMPLEX MODEL 149.7 kg/hr ( 330 Ibm/hr) 


(2) Nodes 288 

(3) Constant Temperatures 2 

(4) Unique Exponents 0 

(5) Temperature Dependent Conductances 0 

(6) Temperature Dependent Heat Inputs o 

(7) Computational Accuracy .0100 

(8) Starting Temperature 25.0 

Are these inputs correct (Y/N) ? Y 
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BRANCH 

NODE 

TON TAG 1. C. S 

BRANCH 

NODE 

TON 

TAG ■ 

I. c, s ' 

403 


249 1 4 

S10 

453 i 


301 

1 

' K224 

404 


301 


K199 


454 1 

118 1 

119 

1 

' S67 

405 


301 


K200 


455 1 


262 

4 

S10 1 

406 

106 

107 1 


S9 1 


456 


301 


K225 

407 


250 


310 


457 


301 


K226 1 

408 


301 1 


K201 


458 1 

119 

120 


S71 1 

409 


301 1 


K202 


459 


263 


S72 

410 

107 

108 


39 


460 


301 


K227 

411 


251 

4 

310 


461 


301 


K228 1 

412 


301 

1 

K203 


462 

120 ' 

121 


S71 

413 


301 

1 

K204 


463 


264 


S72 


108 

109 


39 


464 


301 


K229 

415 


252 


310 


465 


301 


K230 

416 


301 


K205 


466 

““121 1 

122 


S71 

417 


301 


K206 


467 


265 


S72 

418 

109 

110 


39 


468 


301 


K231 

419 


253 


310 


469 


301 


K232 

420 


301 


K207 


470 

122 

123 


S71 

421 


301 


K208 


471 


266 


S72 

422 




39 


472 


301 


K233 

423 


254 


310 


473 

123 

124 


S71 

424 


301 


K209 


474 


267 


S72 

425 


301 


K210 


475 


301 


K234 

426 


112 


39 


476 

124 

125 


S67 

427 


255 


310 


477 


268 

" 4 

S72 

428 


301 

1 

K211 


478 


301 

1 

K235 

429 


301 

1 

K212 


479 

125 

126 

1 

S9 

430 

112 

113 

1 

39 


480 


269 


S10 

431 


256 


310 


481 


301 


S22 

432 


301 

1 

K213 


482 


301 

1 

' K236 

433 


/ 301 


K214 


483 

126 

127 


S9 

434 

113 

114 


39 


484 


270 

4 

S10 

435 


257 

4 

310 


485 


301 


S26 

436 


301 

1 

K215 


486 


301 


K237 

437 


301 

1 

K216 


487 

127 

128 

1 

39 

438 


115 


39 


488 


271 

4 

310 

439 


258 

4 

310 


489 


301 

1 

330 

440 


301 


K217 


490 


301 


K238 

441 


301 


K218 


491 

128 

129 


39 

442 

115 

116 


39 


492 


272 


310 

443 


259 

4 

310 


493 


301 


334 

444 


301 

1 

K219 


494 


301 


K239 

445 


301 


K220 


495 

129 

130 


39 

446 

116 

117 

1 

39 


496 


273 


310 

447 


260 


310 


497 


301 


338 

448 


301 


K221 


496 


301 


K240 

449 


301 


K222 


499 

130 

131 


39 

450 

117 

118 

1 

39 


500 


274 


310 

451 


261 

4 

310 


501 


301 

1 

342 

452 


301 

1 

K223 


502 


301 

1 

K241 


104 
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690 287 286 5 .174E+03 

691 288 287 5 .174E+03 
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TASS GENERAL INPUT MENU - SI Units 

(1) Casa Title: 

TALSR(METRIC)—RUN 5. COMPLEX MODEL, MASS FLOW OP 362.9 kg/hr (300 Ibm/hr) 


(2) Nodas 2S8 

(3) Constant Tanparatiiras 2 

(4) Unique Exponents 0 

(5) Temperature Dependent Conductances 0 

(6) Temperature Dependent Heat Inputs 0 

(7) Computational Accuracy .0100 

(8) starting Temperature 25.0 

Are these inputs correct (Y/N) ? Y 
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.146E+01 

.701E+00 

.240E-01 



,160E+01 

.701E+00 

.910E-01 

.178E+01 

.576E+00 

.178E+01 



576E+00 

174E+03 

174E+03 

174E+03 

174E+03 

174E+03 

174E+03 

174E+03 

174E+03 

174E+03 



.174E+03 

.174E+03 




5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 


5 .174E+03 
5 .174E+03 



5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 



5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 





5 .174E+03 
5 ,174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 



5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
5 .174E+03 
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APPENDIX E. TASS NODAL TEMPERATDSE OUTPUT 


The following output sununarizes the Steady State Thermal 
Analyzer output for the mass flow rates considered in this 
analysis. 
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TALSR(METRIC)- 

-RUN 

1. 

SIMPLE 

MODEL. 

HASS 

FLOW 149.7 

kg/hr 

(330 Ibm/hr) 


Temperatures, 

degc 















1 

25. 

19 

2 

25. 

26 

3 

25. 

40 

4 

25. 

67 

5 

26. 

22 

6 

26. 

,60 

7 

26. 

97 

8 

27. 

31 

9 

27. 

60 

10 

27. 

81 

11 

27. 

97 

12 

28. 

,07 

13 

28. 

13 

14 

28. 

16 

15 

28. 

15 

16 

28. 

10 

17 

28. 

01 

18 

27. 

87 

19 

27. 

65 

20 

27. 

33 

21 

26. 

85 

22 

26. 

58 

23 

26. 

55 

24 

27. 

19 

25 

27. 

56 

26 

27. 

76 

27 

27. 

95 

28 

28. 

03 

29 

28. 

08 

30 

28. 

10 

31 

28. 

12 

32 

28. 

13 

33 

28. 

14 

34 

28. 

15 

35 

28. 

16 

36 

28. 

17 

37 

28. 

18 

38 

28. 

19 

39 

28. 

19 

40 

28. 

19 

41 

28. 

17 

42 

28. 

.12 

43 

27. 

99 

44 

27. 

61 

45 

27. 

15 

46 

27. 

66 

47 

27. 

93 

48 

28. 

05 

49 

28. 

19 

50 

28. 

25 

51 

28. 

29 

52 

28. 

31 

53 

28. 

,32 

54 

28. 

,34 

55 

28. 

35 

56 

28. 

36 

57 

28. 

37 

58 

28. 

37 

59 

28. 

38 

60 

28. 

,39 

61 

28. 

40 

62 

28. 

40 

63 

28. 

41 

64 

28. 

40 

65 

28. 

,37 

66 

28. 

,84 

67 

29. 

28 

68 

29. 

43 

69 

29. 

36 

70 

28. 

87 

71 

28. 

,39 

72 

27. 

,84 

73 

27. 

85 

74 

28. 

42 

75 

28. 

92 

76 

29. 

43 

77 

29. 

,53 

78 

29. 

.39 

79 

28. 

97 

80 

28. 

53 

81 

28. 

58 

82 

28. 

61 

83 

28. 

.62 

84 

28. 

,64 

85 

28. 

,65 

86 

28. 

66 

87 

28. 

,67 

88 

28. 

,67 

89 

28. 

,68 

90 

28. 

,69 

91 

28. 

,70 

92 

28. 

71 

93 

28. 

,71 

94 

28. 

,71 

95 

28. 

,69 

96 

28. 

,64 

97 

28. 

,53 

98 

28. 

43 

99 

28. 

,19 

100 

27. 

,73 

101 

28. 

,18 

102 

28. 

,57 

103 

28. 

.71 

104 

28. 

,77 

105 

28, 

,80 

106 

28. 

,82 

107 

28. 

.84 

108 

28. 

,85 

109 

28. 

. 86 

110 

28. 

,87 

111 

28. 

,88 

112 

28. 

,88 

113 

28. 

.89 

114 

28. 

.90 

115 

28. 

,90 

116 

28. 

,89 

117 

28. 

,87 

118 

28. 

.81 

119 

28. 

,65 

120 

28. 

.48 

121 

28, 

.15 

122 

27. 

,57 

123 

27. 

,61 

124 

27. 

.88 

125 

28. 

,34 

126 

28. 

.64 

127 

28. 

.86 

128 

29. 

,01 

129 

29. 

,11 

130 

29. 

,17 

131 

29, 

.20 

132 

29. 

.19 

133 

29, 

.15 

134 

29. 

,07 

135 

28, 

,95 

136 

28. 

.77 

137 

28, 

.53 

138 

28, 

.23 

139 

27, 

.91 

140 

27. 

,57 

141 

27, 

,09 

142 

26, 

.84 

143 

26, 

.72 

144 

26. 

.66 

145 

25, 

.00 

146 

25. 

,00 

147 

25. 

,00 

148 

25. 

.00 

149 

25, 

.01 

150 

25. 

.02 

151 

25, 

.02 

152 

25, 

,03 

153 

25. 

.04 

154 

25, 

,05 

155 

25, 

.07 

156 

25. 

.08 

157 

25, 

.09 

158 

25, 

,10 

159 

25. 

.12 

160 

25. 

.13 

161 

25, 

.14 

162 

25, 

.15 

163 

25, 

.16 

164 

25. 

,17 

165 

25. 

.18 

166 

25, 

.18 

167 

25, 

.19 

168 

25. 

.20 

169 

25 

.21 

170 

25, 

.23 

171 

25, 

,24 

172 

25. 

.25 

173 

25, 

.26 

174 

25, 

.27 

175 

25 

.28 

176 

25, 

,29 

177 

25. 

.31 

178 

25. 

.32 

179 

25, 

.33 

180 

25, 

.34 

181 

25 

.35 

182 

25, 

.36 

183 

25, 

.37 

184 

25, 

.39 

185 

25 

.40 

186 

25. 

.41 

187 

25 

.42 

188 

25, 

.43 

189 

25, 

.44 

190 

25, 

,45 

191 

25 

.46 

192 

25, 

.48 

193 

25 

.49 

194 

25, 

.50 

195 

25 

,51 

196 

25 

.52 

197 

25 

.53 

198 

25, 

.54 

199 

25 

.56 

200 

25, 

.57 

201 

25 

.58 

202 

25 

.59 

203 

25 

.60 

204 

25, 

.61 

205 

25 

.62 

206 

25 

.63 

207 

25 

.65 

208 

25 

.66 

209 

25 

.67 

210 

25, 

.69 

211 

25 

.70 

212 

25, 

.71 

213 

25 

.73 

214 

25 

,74 

215 

25 

.75 

216 

25, 

.77 

217 

25 

.78 

218 

25 

.79 

219 

25 

.81 

220 

25 

.82 

221 

25 

.83 

222 

25, 

.84 

223 

25 

.86 

224 

25 

.87 

225 

25 

.88 

226 

25 

.90 

227 

25 

.91 

228 

25 

.92 

229 

25 

.93 

230 

25 

.94 

231 

25 

.95 

232 

25 

.96 

233 

25 

.97 

234 

25 

.98 

235 

25 

.99 

236 

26 

.00 

237 

26 

.02 

238 

26 

.03 

239 

26 

.04 

240 

26 

.05 

241 

26 

.06 

242 

26 

.07 

243 

26 

.08 

244 

26 

.09 

245 

26 

.10 

246 

26 

.11 

247 

26 

.12 

248 

26 

.14 

249 

26 

.15 

250 

26 

.16 

251 

26 

.17 

252 

26 

.18 

253 

26 

.19 

254 

26 

.20 

255 

26 

.21 

256 

■ 26 

.22 

257 

26 

.23 

258 

26 

.24 

259 

26 

.25 

260 

26 

.26 

261 

26 

.27 

262 

■ 26 

.28 

263 

26 

.30 

264 

26 

.31 

265 

26 

.32 

266 

26 

.32 

267 

26 

.33 

268 

26 

.34 

269 

26 

.35 

270 

26 

.36 

271 

26 

.37 

272 

26 

.38 

273 

26 

.39 

274 

26 

.40 

275 

26 

.41 

276 

26 

.42 

277 

26 

.43 

278 

26 

.44 

279 

26 

.45 

280 

26 

.46 

281 

26 

.47 

282 

26 

.48 

283 

26 

.48 

284 

26 

.49 

285 

26 

.49 

286 

26 

.49 

287 

26 

.49 

288 

26 

.49 

301 

40 

.00 

302 

25 

.00 
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